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Abstract. We discuss the application of the momentum-sheU renormahzation group 
' method to the interacting homogeneous Bose gas in the symmetric and in the 

symmetry-broken phases. It is demonstrated that recently discussed discrepancies 
are artifacts of not taking proper care of infrared divergencies appearing at finite 
^ ' temperature. If these divergencies are taken into account and treated properly by 

, means of the e-expansion, the resulting renormahzation group equations and the 

' corresponding universal properties are identical in the symmetric and the symmetry- 

I broken phases. 
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PACS numbers: 03.75.Hh,05.30.Jp,64.60.Ak 
1. Introduction 



O . Renormalization techniques have been employed for the study of interacting Bose gases 
^ ' near the critical temperature, because in this temperature regime the fluctuations 
dominate the mean field. The bulk of this work was written before the experimental 
^ ! realization of Bose-Einstein condensation (BEC) in ultracold atomic gases [11121131110 
IHIIZIEI- It was shown that only at zero temperature the quantum nature of the three- 
dimensional Bose gas differentiates it from a three-dimensional two-component classical 
field theory. At any finite temperature the Bose gas converges to the classical theory as 
the fixed point of the renormalization group (RG) equations is approached. 

Therefore the calculation of the universal properties of the Bose gas can be 
performed in the same way as for a classical theory, that is through the e-expansion 
where e = 4 — D and D is the number of spatial dimensions, e.g. jHllTOlEll^- In 
D = 3, although the results of the expansion up to second order in e are in remarkable 
agreement with experimental values of critical exponents (measured in He^ experiments, 
but due to universality applicable in the case of Bose gases as well [IH|), higher-order 
results diverge from the experimental values. The reason is that the £-expansion is 
asymptotic, as first noted in and to obtain meaningful results when higher orders 
in e are included, one has to make use of resummation techniques, see e.g. [T^. This way 
critical exponents have been calculated up to fifth order in [Tni[IZl[IHl[inil2nil2Ill22| , see 
also for corrections j2SlEl] and improvements [2HJ- However, the results thus obtained 
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are sensitive to the way the resummation is performed and consequently somewhat 
ambiguous, see e.g. P5] . 

An alternative method is to calculate the universal properties perturbatively as 
series of powers of g* {g* being the infrared stable fixed point for the interaction g) 
directly in D = 3, as first suggested in These series are then truncated to order 
g*^ where L is the number of loops in which the calculation is performed. Though this 
method is fundamentally less satisfactory than the e-expansion, see e.g. ^3], it can be 
used in the regime of small but non-zero chemical potential. It has been employed for 
the calculation of critical exponents up to seventh order in g* for = 0, 1, 2, 3 pH|l29p26j 
and for arbitrary jHIUnT] where A^ is the number of components of the vector field. 
The series in g* are again asymptotic and have to be resummed. There is in general 
agreement with the corresponding e-expansion results. We will be referring to this 
technique as the direct method. 

After the experimental realization of BEC in ultracold atomic gases 
because of the renewed interest in these systems, a new generation of papers on the 
renormalization of Bose gases appeared. Starting with |35], a series of papers relied 
on the so-called momentum-shell approach p ^ l!?7 | l!?H | l39j . In this method, momentum 
shells around the cutoff are successively integrated out directly at D = 3 according to 
Wilson's method, but unlike in the direct method no expansion of the critical exponents 
over g* is performed. 

This apparently new method, when applied in the symmetric (normal) phase, yields 
universal results which, when compared to experimental values, are worse than even the 
first-order e-expansion results. However, when the momentum-shell method is applied 
to the symmetry-broken phase, it yields results which are far better than the first- 
order e-expansion and, in fact, as good as the results of the second-order e-expansion. 
Based on this observation, it was assumed that the reliability of the momentum-shell 
method increases when it is used in the symmetry-broken phase, and a calculation of 
non-universal properties (for example transition temperature versus scattering length) 
from the symmetry-broken phase RG equations was attempted. For this reason, one 
may now wonder if applying the e-expansion or the direct method in the symmetry- 
broken rather than in the symmetric phase as is usually done would improve the results 
of these methods. 

Here we show that the distinction introduced by the momentum-shell method 
of ISSllSSllSZllSHllSni between RG methods applied in the symmetric and the symmetry- 
broken phases is an artifact of not taking care of the infrared divergence which appears 
at finite temperature in the interacting Bose gas theory. When this divergence is taken 
into account and treated properly, as in a classical theory, by means of the e-expansion, 
the RG equations and the resulting universal properties are identical in the symmetric 
and in the symmetry-broken phases. 

Furthermore, even if one focuses on the regime of small but non-zero chemical 
potential, where the direct method applies, the use of the direct method in the 
symmetry-broken rather than the symmetric phase deteriorates the results instead of 
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improving them. 



2. Derivation of RG equations 

The partition function of the homogeneous s-wave interacting Bose gas is 

/5, V, g) = Tre-^i^-^'^) = J 0*]e-^['^''^*l (1) 

with the Euchdean action 

1 r'^P r ( 



dr 2m ^ 



(r,x) + ||</>(r,x)|4. (2) 



We give an outhne of the basic steps of the renormahzation procedure. More details 
can be found in, e.g., jSHlEZlEH]- 

In order to implement the first step of the RG procedure (Kadanoff transformation), 
we split the field (f){x) into a long- wavelength component 0<(x) (slow field) and a short- 
wavelength component 6(j)y{x) (fast field). The fast field involves Fourier components 
which are contained only in an infinitesimally thin shell in momentum space of thickness 
Ae~' < IpI < a near the momentum cutoff A, whereas the slow field has all its Fourier 
components in the sphere whose center is at the origin of the momentum space and 
whose radius is Ae~'. 

We now perform the one-loop calculation of the effective theory of the slow field. 
We integrate out the fast field and expand the resulting effective action in powers of g 
keeping up to order g"^. This perturbative effective action is equal to the original action 
(0) (with the field replaced by the slow field 0<) plus two additional terms. The first of 
these additional term is proportional to g and therefore quadratic in the modulus of the 
slow field. The second additional term is proportional to g^ and quartic to the modulus 
of the slow field. The effective Lagrangian of this theory can be cast in the form of 
the original Lagrangian because the additional terms produced by the integration over 
the fast fields have such a form that they can be considered as corrections to fi and g. 
Thus, after one infinitesimal integration, the chemical potential and the interaction in 
the effective action are 

li' = 11 + dii = 11- gdo / dp E{p) - fi], 

J &Vp 

g' =g + dg = g- g^dn f dp p''-'f2[P, E{p) - /i] (3) 

where E{p) = p^ /2m and 

/i(/3,E) = 2Ar(/3E) + l, 

/2(/3, E) = ^ + ^^^^^^^ + ApNiPE) [1 + N{PE)] (4) 

with the Bose-Einstein distribution N{x) = l/(e^ — 1). The density of states do can be 
expressed in terms of the surface Qo of a unit hypersphere in D dimensions according 
to 

1 27r^/2 
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The infinitesimal momentum shell which is integrated out is denoted by 6Vp. We note 
that the integration procedure has no effect on the inverse temperature or the slow field, 
i.e., 

P' = = </)<. (6) 

Now the only remaining difference between the form of the effective action and that of 
the original action is that in the original action the momentum of the field is integrated 
from to the momentum cutoff A whereas in the effective action the momentum of the 
slow field |q'| is integrated from to Ae~K This difference is eliminated by a trivial 
rescaling of the momentum \q{l)\ = |q'|e' which induces a trivial rescaling of the rest of 
the parameters of the effective action, 

^(/) = /i'e^', gil) = g'e^'-^^\ = /3e-'\ 0<(/) = 0<e™/l (7) 

In terms of the rescaled parameters, equations assume the form 

= e^V - e^'gdo T dp p^'Vil/^, E{p) - /i], 

^7(0 = e(2-^)'^ - e(2-^Vdz. r dpp''-'MP,Eip)-iA. (8) 

The system of (jHl) becomes autonomous if one solves it together with 

= e-''p. (9) 

We perform the integrations in (jHJ over the infinitesimal momentum shell SVp and keep 
terms only up to first order in /, i.e., 

fi{l) = /i + 2^/ - gdoA^fi [(3, E{A) - fi]l, 

gil) =g + i2- D)gl - g'dnA''f2[P, i?(A) - /i]^ (10) 

Repeating the above procedure of integrating out shells of high-momentum and 
rescaling, we find the RG equations for the chemical potential and the interaction, 

^ = 2/i(/) - gil)dnA''fimi),E{A) - 

^ = (2 - D)g{l) - g{lfdj,A^h[m.E{A) - m (H) 

where the number of renormalization steps / is a continuous parameter running from 
(no shells integrated out) to = oo (all shells integrated out). The trivially rescaled 
quantities defined in (jZj) also become functions of I. 

Setting (7 = and consequently g{l) =0 in (fTT|) one can study the ideal Bose gas 
both at zero (/3 = oo) and at finite temperature (/3 ^ oo), e.g. jHEI- One can also study 
the zero-temperature, interacting Bose gas, e.g. jlHlllIlElIll] • Finally, one can study the 
problem in its full generality, the finite-temperature interacting Bose gas. 

In reference jH], the symmetric phase of the finite-temperature interacting Bose 
gas in an arbitrary number of dimensions was considered. It was shown that, for any 
finite temperature, the interacting Bose gas is driven towards a two-component classical 
system, as the fixed point of the RG equations is approached, see also 01211311100111 • 
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Our formalism is different from that of 8J or 1 , ^HSHHElinilZj , so it may be of some 
interest to show that we can come to the same conclusion. In the course of the following 
calculation we will also see how the finite temperature theory we are examining here 
develops an infrared divergence. 

The RG method which we have been using is perturbative over the interaction g. 
Therefore, as in the zero-temperature interacting gas case, e.g. we examine the 
behaviour of (fTT|) near the fixed point = (0,0) of the unperturbed system, that 

is the finite-temperature ideal gas. Near this fixed point the Bose-Einstein distribution 
can be expanded as follows 

and therefore equations (jH)) near the fixed point become 

f-A 1 



/x(/) = e^V - (^^^gdo I dp p^-^ 2 

JAe-' 



We observe that, at (/i = 0,/? = 0), the last term in each of the above equations is 
infrared divergent for all p and all D. This divergence can be treated by redefining 
variables; we recast the above equations in terms of the new variable g = g/P, 

fi{l) = e^V - e^'~gf3dD ^ dp p^'^ f\[(3, E{p) - /x], 

m = e^'-^^'-g - e^'-''^'f(3do l"^ dp p^"^ h[(3. E{p) - /x] (14) 

JAe-' 

which can be also written in a differential form 

2/i(/) - mmdoP^'^hm). m) - 



di 

^9(1) _ (4 _ jjyg^i^ _ ~g{if(3{l)dj,K^U[(3{l), E{A) - /.(/)]. (15) 



dl 

These are exactly the equations (4.6) of |H] for the chemical potential and the interaction 
supplemented by equations (4.12) of jS] in the classical regime. 

Near the unperturbed fixed point (/i*,/9*) = (0,0), Eqs. (fT^ reduce to 



cA 1 



= e^V - e^^gdo dp p^ _ 

JAe-' E[p) — fl 

Ml) ^ e<'-.', - .<--^,=.. * 5^^. (16, 

We note that apart from some numerical coefficients, equations fjl6|) are identical 
with the equations ensuing from the classical Landau- Ginzburg- Wilson (LGW) reduced 
Lagrangian, see, e.g., equations (6.40)-(6.42) of [15] . 

The different coefficients are due to the fact that the Bose gas theory involves a 
complex field, whereas the classical LGW theory uses a two-dimensional vector field. As 
we will see later, this difference is insignificant in the sense that the two theories have 
the same universal behaviour. 
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Finally, it is convenient to cast (jl5|) in terms of dimensionless variables M = f3\iJ,, 
G = m'^A^~'^gdr)/h'^, = P\E{A) = 1/2 (so that we keep track of the energy terms), 
and b = (3/ /3a, with /?a = m/{%^!s?'), which yields 
dM{l) 



dl 



2M{l)-Gmi)fi[m,Ey-M{l)], 



^ = (4 - D)G{1) - G{lfb{l)h[b{l), E^ - M{1)]. (17) 

It may not be immediately obvious that the variable G is dimensionless. However 
recalling that g = VLoh^a^~'^ /m, see jH], we can rewrite it as G = ^7/)(aA)'^^^/6 which 
is clearly dimensionless. 

2.1. e-Expansion 

Since equations (fTBj) have the same structure as the classical RG equations, the same 
difficulties in the infrared regime appear. In particular around /x = and for D < A the 
integral in the equation for the interaction is infrared divergent. It is interesting to note 
that this divergence originates in the terms 

E{p) - /i 

of f2[P, E{p) — fj], see equation (jl)). As in the classical case, this divergence is cured by 
performing the e-expansion. We identify 4 — D in the equation for the interaction with 
e. Furthermore we assume that fi and g (and therefore M and G) are of the same order 
and expand the RG equations (jl7|) up to second order in these variables which yields 

dM{l) 



dl 



2M{1) - G{l)b{l) {2N[b{l)Ey] + 2b{l)N[b{l)E^]{l + N[b{l)Ey])M{l)} 



^ = eG{l) - G{lfb{l) I ^ + 2Aa&(^)^>] + Ab{l)N[b{l)E,]{l + NmE>])] ■ (18) 

This system has a trivial fixed point at (M*, G*) = (0, 0) with eigenvalues 
Ai = 2, 

\^=£ = 4:-D. (19) 

Therefore, the eigenspace of A2 corresponds to the unstable direction for D < 4, to the 
marginal one for D = 4, and to the stable direction for D > 4. There is also a non-trivial 
fixed point 



(M*,G* 



5e 



10 -2e' (10 + 2e)2 
Up to first order in e, its eigenvalues are 

Ai = 2 - le, 
5 

A2= - e = D-4, (20) 

and consequently the eigenspace of A2 corresponds to the stable direction for D < 4, to 
the marginal one for D = 4, and to the unstable direction for D > 4. 
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For the case of physical interest, D = 3, we set e = 1 in the expressions which 
we have aheady expanded up to first order in e. A simple example of a universal 
property that we can now calculate is the critical exponent for the correlation length 
z/ = 1/Ai = 0.600 + Oi^e^). This is exactly the same as the result found in the e- 
expansion study of a classical two-component LGW theory in D = 3, e.g. piTIHITT]. see 
also |7j. A finite temperature interacting Bose gas in three spatial dimensions belongs 
to the same universality class as a two-component classical field theory in three spatial 
dimensions (Zj. 



2.2. Direct Method 

We set D = 3 directly in (fTTj) or equivalently in (fTTj) . The fixed point is 

{^^\~g*) = {E{^)/QM'E{^f/{l%^^)) 

or equivalently in dimensionless units {M*,G*) = (1/12,5/144). We then linearize 
around the fixed point and calculate the eigenvalues 



Ai = — -1728G* + V2985984G'*2 + 57600^* + 625 + 75 
50 

1 



A2 = — -1728^* - V2985984G*2 + 576OOG* + 625 + 75 . (21) 
50 V / 

Finally, we expand the critical exponent z/ = 1/Ai up to first order in G* because we 

are performing a one-loop calculation, 

v = - + + 0(G*2) = 0.600 + 0(G*2). 

This is the same result as for a two-component classical theory, see e.g. PUI I^. 



2.3. Momentum-Shell Method 

This approach is used in [22], see also [SZUSHl, and corresponds to setting D = 'i directly 
in ()11|) or equivalently in ()17p as in the direct method. In this case the RG equations 
(|11|) coincide with equations (11a) and (lib) of The fixed point is the same as in 
the direct method. However now we set G* = 5/144 in the eigenvalues 
_ 3 + V249 
10 ' 

3 - 7249 , , 

A2 = 22 

10 ^ ^ 

and without expanding we find u = 1/Xi ~ 0.532. 



3. Symmetry-broken phase 



If we spontaneously break the global U(l) symmetry of (j21) by introducing the most 
probable configuration 0, we find 



-pv 



nig 
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h 



d h 



/i----— V^-/i + 2gno 
OT 2m 



•x) 



/■ftp r 

/ dr / rf^x0*(r,x) 

JO Jv 

/ dr / f/^x [0*(r, x)0*(r, x) + 0(r, x)0(r, x)] 

JO Jv 

/ rfr / f/^x [0*(r, x)0*(r, x)0(r, x) + 0*(r, x)0(r, x)0(r, x)] 

JO Jv 



+ |- dr c?^x 0*(r, x)0*(r, x)0(r, x)0(r, x) (23) 

where no = |0p = fi/g is the condensate density, and 0(r, x) now denotes the fluctuation 
around the most probable configuration (f). 

It is possible to write down the RG equations for the symmetry-broken phase. The 
calculation is significantly more complicated than in the symmetric phase, for details on 
the derivation of the RG equations see 13^1 for D = 3 and [27|EH] for arbitrary D. We 
focus here on the RG equations for the chemical potential and the interaction 

. um - c(OMO {?^M)M)!|2iv(K0A(0) + 11-1 

+ ^^''''^i-,;"^''''^ (Oiv(KOA(0)|iv(MOA(0) + i|}. 

^ = (4 - D)G(l) - Oilfbil) |i?^=ifp!|2/V(M()A(()) + 1] 

+ '^\^y''"' ''(')^(K')A(i))|/V(M')A(0) + 1]} (24) 



where A(/) = y-El + 2M(/)i?>. We note that the above RG equations coincide with 
(jl7|) of the symmetric phase for M = 0. 

3.1. e-Expansion 

At the fixed point (yU*,/5*) = (0,0) of the unperturbed system an infrared divergent 
term appears in the course of the derivation of the RG equation for the interaction. 
This term is exactly the same as in the symmetric case, i.e., 

but now originates in the terms 

m^iV(;3*(p,) and ^2EM^^M 



with 6{p) = ^E{pf + 2^iE{p). 

We can easily apply the e-expansion technique to equations in order to cure 
the infrared divergence. As we have seen in the symmetric phase, all we have to do is 
identify A — D with e in the equation for the interaction and then expand the chemical 
potential equation up to first order in M and the equation for the interaction up to 
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zeroth order in M. We thus find the e-expansion RG equations 
dM{l) 



dl 



2M{1) - G{l)h{l) {2N[h{l)Ey\ + 2h{l)N[h{l)E>\{l + N[h{l)E^])M{l)} 



^ = eG{l) - ami) 1 1 + 2AaK0i^>] + 4K/)iV[6(/)i5>](l + N[bil)E,])^ . 

(25) 

Comparing (j^^ to ()18|) we see that they are exactly the same! In other words, when we 
perform the momentum- shell integrations together with the £-expansion, the symmetric 
and symmetry-broken phases yield exactly the same RG equations, and consequently 
identical universal properties. 



3.2. Direct Method 

As in the symmetric phase, it consists of setting directly D = 3 in (jUj). The non-trivial 
fixed point is (M*,(5*) = (1/2,1/4). The eigenvalues of (j21I), when it is linearized 
around the fixed point, are 

Qtt^ - 12G* + V3(-16G*2 + 2A7c^G* + Sn^y^^ 



Ai 



2 



- 12G* - V3{-16G*^ + 2471^0* + 37i^y/^ , , 

= ^2 ' (26) 



and therefore 

zy = 0.500 + 0(G*2). 



3.3. Momentum- Shell Method 

We set directly D = 3 in as in the direct method. The resulting RG equations, 
when recast in the dimensionful variables fi and g, coincide with (29a), (29b) of jHS], see 
also equations (16), (17) of [37J. We now substitute G* in the eigenvalues of the direct 
method 

3 + ^33 



Ai 



6 



A2 = (27) 
and without expanding we find u = 1/Ai ~ 0.686. 



4. Comparison and Conclusion 

The first-order e-expansion and the momentum-shell method used in and in 
subsequent papers [36jl37p38j do not yield identical universal properties. As we have seen, 
the first-order e-expansion gives the value u = 0.600 + 0{e'^) for the critical exponent 
of the correlation length both in the symmetric and in the symmetry-broken phases. 
However the momentum-shell method gives u = 0.532 in the symmetric and u = 0.686 in 
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the symmetry-broken phase. These results are to be compared to the experimental value 
V = 0.670, see e.g. fT^. It seems as if the momentum-shell method in the symmetry- 
broken phase yields the best result (in fact, almost as good as the second order in the 
e-expansion result u = 0.655 + 0{e^), see e.g. |15|). However, the results do not always 
justify the method used to obtain them. For example, in subsection 2.1, had we not first 
expanded z/ = 1/Ai up to first order in e and then set e = 1 but used the unexpanded 
formula (j^UI) for Ai, we would have found u = 0.625. Although this is closer than 
u = 0.532 to the experimental result u = 0.670, it is clearly incorrect. 

The main point of this paper is that the momentum- shell method results depend 
on whether the RG calculation is performed in the symmetric or the symmetry- 
broken phase, whereas the first order e-expansion results do not. The dependence of 
the momentum-shell results on the phase is an artifact of not avoiding the infrared 
divergence of the Bose gas theory. Because for any finite temperature the Bose gas 
theory has the same infrared behaviour as a classical two-component theory, the e- 
expansion as we know it from classical papers, e.g. [TTj, can cure the infrared divergence 
and yield reliable results. 

Furthermore we note that, even if we do not worry about the infrared divergence, 
and use the direct method, applying it in the symmetry-broken phase makes the result 
for the critical exponent of the correlation length worse rather than improving it as one 
would have expected from the momentum-shell method. 

The above discussion and comparison of values does not take into account the effect 
of several approximations (derivative expansion, polynomial truncation of the effective 
action, assumptions about the sharpness or smoothness of the infrared cutoff separating 
the fast from the slow modes) which were employed in the course of the calculation. For 
discussions on these approximations see |i31l361l39j . However since these approximations 
were used in the part of the calculation which is common for all three methods, that is 
the derivation of the RG equations (fT7j) (or equivalently (fTTjl ) and (j21I), a comparison 
between these methods is still valid. 
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